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Abstract 

We study the diffusive transport of neutrinos in a newly born neutron 
star to explore its sensitivity to dense matter properties. Energy and 
lepton number which are trapped during the catastrophic implosion 
diffuse out on the time scale of a few tens of seconds. Results for different 
dense matter models are presented. 

1 Introduction 

The core of a massive star implodes when its mass exceeds the Chandrashekar 
mass. The hot and dense remnant formed subsequent to the implosion is a 
protoneutron star. Numerical simulations of the implosion (and the subsequent 
formation of a shock wave at core bounce) indicate that, due to the high 
densities and temperatures, most of the star's gravitational binding energy 
and lepton number released remains trapped in the star as neutrinos. The 
general features of the early evolution have been discussed in prior work |], § • 
The object of this work is to elucidate the role played by the microphysical 
inputs (equation of state (EOS) and neutrino opacities) on the macrophysical 
evolution of the protoneutron star. 

2 The Basic Equations 

The structure of the star is assumed to be in quasi-static equilibrium, and 
is described by the general relativistic equation for hydrostatic equilibrium. 
The loss of lepton number and energy due to neutrino flows is treated in the 
diffusion approximation. The equations governing the early evolution may be 
found in [|I|, |2|. To solve the structure equations we need to specify the finite 
temperature EOS discussed in §3. Other microphysical inputs required to solve 
the transport equations are the neutrino mean free paths and bulk properties 
such as the specific heat and the nuclear symmetry energy. 



2.1 Deleptonization 

To illustrate how microphysical inputs influence lepton number transport, we 
examine the Newtonian diffusion equation || 
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Here Yl = Y e + Y v is the total lepton fraction and \x v is the local neutrino 
chemical potential. The neutrino mean free path \{E V ) enters Eq. (1) via the 
unnormalized diffusion coefficients D n (and D n for anti-neutrinos) defined by 
D n = Jo°° dE v \ V (E V ) f v (E v )(l — f u (E u )). The time scales crucially depend 
on the magnitude of \(E U ) and its energy dependence. For electron neutrinos, 
the dominant contributions to X(E U ) arise from the charged current reaction 
v + n — > e~ + p. The rate of this reaction is sensitive to both temperature and 
the composition of the ambient matter [3-6]. Under degenerate conditions, the 
phase space for this reaction is proportional to fi n T 2 (£t + Ej), where \i n is the 
neutron chemical potential and fi = fi n — fi p . The left hand side of Eq. (1) 
illustrates one way the EOS directly influences the lepton number flows. Note 
that = fp^^Sf , and that j& is directly related to the nuclear symmetry 
energy S(n B ). Since fi = 4S(rib)(l — 2Yp), where Y p is the proton fraction, the 
nuclear symmetry energy will play a crucial role in the deleptonization phase. 
Further, the possible presence of hyperons at high density affects both |^ and 
the neutrino mean free path A will be altered 0, || . 

2.2 Cooling 

At early times, the flow of electron neutrinos due to chemical potential gradi- 
ents results in net heating rather than cooling. This Joule heating and negative 
temperature gradients warm the central regions until the center reaches the 
maximum temperature. At this stage, the center begins to cool as the thermal 
gradients force all six neutrino species radially outwards. Due to the high den- 
sities and temperatures present, the pair production processes are sufficiently 
rapid to ensure that /i and r neutrinos are in thermal equilibrium. Energy 
transport is governed by 
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All general relativistic corrections and terms arising due to lepton number 
gradients have been dropped to highlight the essential role of the microphysics. 
The dominant opacity for // and r pairs is that due to scaterring on baryons and 



electrons, as the charged current reactions are kinematically inaccessible. The 
charged current reactions are still the dominant source of opacity for electron 
neutrinos; thus, relative to the \x and r pairs their contribution to energy 
transport is significantly suppressed. Under partially degenerate conditions, 
the phase space for scattering reactions is roughly proportional to M* 2 T 3 (note 
that the ambient temperature determines the average E u ); M* arises since it 
determines the density of states at the Fermi surface. The left hand side 
of Eq. (2) is directly related to the specific heat since Tds = C v dT, and 
is strongly dependent on the effective baryon mass M*. Neutrino mean free 
paths depend on the same physics, since they are proportional to the number 
of target particles which do not suffer significant Pauli blocking. 

3 Microphysics 

3.1 The Equation of State (EOS) 

The dense matter EOS is an important ingredient in protoneutron star sim- 
ulations. For a detailed discussion of a variety of dense matter models see 
[0]. In this work, we have explored four different models (see Table) to study 
the early evolution. We broadly classify them as soft (GM3) and stiff (GM1) 
models ||, and models can contain (GM1-H, GM3-H) or not contain (GM1-N, 
GM3-N) hyperons. Hyperons begin to appear only towards the late stages of 
the deleptonization 0. This results in these models having significantly lower 
maximum masses for Y u=0 (4 th column) and the existence of a range (3 rd & 4 th 
columns) of initial masses that become unstable upon deleptonization. 
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M^tUT = 0) 


M-^(T = 0) 


GM1-N 


n,p + leptons 


2.52 


2.76 


GM1-H 


n,p,H + leptons 


2.25 


1.96 


GM3-N 


n,p + leptons 


2.16 


2.35 


GM3-H 


n,p,H + leptons 


1.95 


1.74 



3.2 Neutrino Interactions 

Neutrino interactions are significantly modified in a hot and dense medium. 
The effects of Pauli blocking and strong interactions on the weak interaction 
rates have been investigated in Refs. [3-6] which showed that density and tem- 
perature dependence of the neutrino mean free path are significantly different 
from those employed in earlier proto- neutron star simulations [l], [2|. Strong 
interaction effects may be incorporated by employing appropriate dispersion 
relations and effects of Pauli blocking may be incorporated exactly. In ad- 
dition, Ref. || shows that the neutrino mean free path also depends on the 
composition of the ambient matter. In particular, the appearance of hyperons 
leads to significant reduction in A„ during the cooling phase. 



4 Results 
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Figure 1: Central temperature and neutrino chemical potential vs time (model GM3). 
Right (left) panels: matter with (without) hyperons. Curves are labelled by the initial 
baryon masses. Models Mb > My v =Q — 1.74Af Q are unstable and collapse as they delep- 
tonize (denoted by an asterix). Note, however, that if the baryon mass does not significantly 
exceed this value the meta-stable state lasts for a very long time. 



We turn now to results of simulations which include the full effects of 
general relativity. In Fig. 1, the time evolution of the central temperatures are 
shown for models with (right top) and without (left top) hyperons for different 
initial masses. The deleptonization times are related to the time evolution 
of the neutrino chemical potential, which are shown for models with (right 
bottom) and without (left bottom) hyperons. The identification of time scales 
associated with deleptonization and cooling depend on their precise definition. 
Here, we opt to define the deleptonization time tjj as the time it takes for the 
central value of \i Ve to drop below 10 MeV, and the cooling time tq as the time it 
takes for the central temperature T c to drop below 5 MeV. Fig. 2 clearly shows 
the generic trends; (1) Softening leads to higher central densities, and thus 
higher temperatures, both of which act to decrease \ u , (2) Hyperons decrease 
the baryon degeneracy in the central regions, and since the neutrinos couple 
quite strongly to the hyperons J4[, X u is reduced; and (3) Hyperonization is 
accompanied by compressional heating and neutrino production, both of which 
delay the cooling and deleptonization times. Hyperons thus always act to 



increase the deleptonization and cooling times. Fig. 2 also shows that a softer 
EOS favors longer diffusion times due to higher temperatures and densities 
in the inner regions of the star. To discriminate observationally between these 
different dense matter scenarios, the neutrino luminosities must be folded with 
the response of terrestrial detectors. 
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Figure 2: Deleptonization and cooling times. Left panel: td for various models as a 
function of the initial baryon mass. Right Panel: tc for the different models. 
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